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Abstract 



The energy of a test particle orbiting a Schwarzschild black hole is 
quantized owing to the quantization of the angular momentum. For small- 
1 est stable circular orbit, the excitation energy is found to resemble closely 

the expression for the temperature of the Hawking radiation. This result 
,1 ^ ■ is consistent with the Unruh effect for orbiting test particle. The predicted 

j^J} energy quantization might be observable by studies of the red-shifted 21- 

cm line of neutral hydrogen orbiting a primordial black hole with mass of 
the order of that of Earth. 

PACS number(s): 04.70.Bw, 04.70.Dy, 95.30.Sf 



1 Introduction 

Long ago, Dirac showed that the existence of magnetic monopoles implies quan- 
tization of electric charge [1]. In 1985, Zee [2] proposed a gravitational analog 
of Dirac quantization condition. This analog is based on the theory of gravito- 
electromagnetism (GEM) [3]. In ordinary GEM, which is the analog of Maxwell 
theory it is the mass of the test particle and the gravitational acceleration that 
play the role of the electric charge and the electric field, respectively. The GEM 
magnetic field is given by the GEM analog of the Biot-Savart law in which the 
source is the matter current density. This field is a divergenceless quantity ev- 
erywhere. Following Dirac, Zee [2] considers the possibility that a GEM analog 
of the magnetic monopole, a "gravitipole" , exists. Consequently, the magnetic 
field of the gavitipole satisfies the equation 

V.B = j5 3 (x) (1) 

Using the Gauss law, we can integrate Eq. (1) to observe that the constant 
7 is equal to the total magnetic flux $ emanating from the gravitipole. The 
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quantization of the mass of the text particle is obtained by expanding the action 
S of the test particle in the gravitational field. Zee [2] focuses on the term 
representing the interaction of the test particle with the vector potential C of 
the gravitipole 



where m is the rest mass of the test particle. As the test particle is moved 
along a closed loop around the gravitipole, its wave function acquires a phase 
Sint/h — —m'f/h. Single-valuedness implies that this phase be equal to 27rn 
where n is an integer. This yields the quantization condition [2] 



Owing to the hypothetical nature of the gravitipole, the constant 7 is not known. 
Only the upper bound of the energy splitting can be estimated from the absence 
of the quantization effects on energy level splittings in atoms and molecules [2] . 

The present effort was inspired by the paper of Zee [2]. We consider a 
quantization of the energy that is due to the gravitational field of a black hole. In 
distinction from Ref. [2] , it is not the rest mass of the test particle that becomes 
quantized. Rather, it is the radius of the particle orbit that is quantized as a 
result of quantizing the angular momentum. For a test particle, moving along a 
circular path around the black hole, the quantization is due to the fact that the 
azimuthal angle is a periodic function of time. Since the canonical momentum 
that is conjugate to this angle is the angular momentum, the corresponding 
action variable is the angular momentum times 2tt. Applying the postulate of 
the old quantum theory we set the action variable equal to nh where n is an 
integer and h is the Planck constant [4,5]. In this way, the energy quantization is 
linked to the quantization of the angular momentum of the test particle orbiting 
a black hole. 

In Sec. 2, we consider a Schwarzschild black hole and derive the Lagrangian 
and the canonical momenta for the test particle. Using these momenta, we 
obtain in Sec. 3 the effective potential and the radial equation of motion of the 
time-like geodesies. From the effective potential, we determine the condition 
that the angular momentum must satisfy in order that the geodesies are stable 
circular orbits. In Sec. 4, we derive the quantization condition for the energy of 
the test particle. We find that the mass quantum is a function of the mass of the 
black hole and the radius of the orbit. For the smallest stable orbit, the formula 
for the excitation energy resembles closely the expression for the temperature 
of the Hawking thermal radiation [6]. In Sec. 4, we attempt to clarify this 
similarity by comparing the period of the orbital motion with the period of the 
Euclidean time. Sec. 5 is devoted to estimates of the energy quantization in 
various astronomical settings. The relation of the excitation energy to Hawking 
temperature is discussed in the Appendix in terms of the Unruh effect [22] . 
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2 Lagrangian and canonical momenta of test par- 
ticle 



We consider a test particle of mass m orbiting a Scharzschild black hole. We 
confine ourselves to circular geodesies in the equatorial plane. Then the metric 
takes the form [7,8] 

ds 2 = -(1 - T -^)c 2 dt 2 + (1 - T -^)- l dr 2 + r 2 d<t> 2 (4) 
r r 

where r$ = 2MG/c 2 is the Schwarzschild radius for a black hole of mass M. We 
assume that M is much larger than the Planck mass, ensuring that rs is much 
larger than the Compton radius, rc — fi/Mc, so that quantum fluctuations of 
the black hole (zitterbewegung) can be disregarded [9]. 

The Lagrangian of the test particle can be expressed in terms of the metric 
components gij as follows [7,8] 

£ = --^9ijX l x 3 = ~ (gooc 2 i 2 + g rr r 2 + g^ 2 ) (5) 
where we adopt the notation 

.n • dt . _ dr i dd> 

x° = ct = c— ,x 7 = — ,i = -r- 6 
dr dr dr 

Since we are considering time-like geodesies, we choose r to be the proper time. 

The canonical momenta, p a = d£/dx a , are obtained with the use of Eqs. 
(4) and (5) yielding 

d £ v 
Pt = — T = ~mc 2 g oa i = mc 2 (l-—)i = E (7) 
at r 

d£ 

P4> = = -mg^4> = -mr 2 4> = -L (8) 
ocp 

Pr = -gT = -mg rr r = -m(l - — ) r (9) 
where E and L are the total energy and angular momentum of the test particle. 



3 Effective potential for radial motion 

The equation for radial motion of a geodesic is obtained from the relation 

g m pl + g rr pl + g^ P l + m 2 c 2 = Q (10) 

where the first three terms describe the magnitude of the energy-momentum 
four-vector (p = E/c,p). Introducing into Eq. (10) the quantities g y from Eq. 
(4), and the canonical momenta from Eqs. (7-9), we obtain [8] 
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E 2 



±L (i_^)-i +m 2 ( i_- 



r 2 

-1 -2 , - L/ 2 2 

V z + — + m z c z 







(11) 



C r r 

Following the notation of Rcf. [7] , we introduce the energy and momentum 
of the test particle per unit rest mass 



m m 

Introducing these definitions into Eq. (11), we obtain 



(12) 



2 2 
TO C 



^r(i-^r 1 + 3(i-^r 1 + 



+ i 



= o 



From this equation, we obtain the equation for radial motion in the form 
a £ 2 1 12 



_-_ 4- 2\ fi _ T _S\ 

2 2c 2 2 V r 2 A r ) 



(13) 



(14) 



Letting G = c = 1, this equation takes the form of Eq. (6.3.14) of Ref. [8]. The 
effective potential for the radial motion is obtained from Eq. (14) in the form 



T / W 

The extrema of this potential are found from the equation 



(15) 



dv e 



eff 



dr 2r 4 
The roots of this equation are 



R± 



L 2 

2 

c rs 



± 



2L 2 

c 2 rs 



L 2 

c z rs 



r + 



3L- \ 



r 3. 



(16) 



(17) 



We see that the roots are real only if L 2 > 3c 2 r|. In this case, R + is a true 
minimum of the effective potential. This can be verified by computing, using 
Eq. (15), the second derivative of the effective potential which turns out positive 
at r = R+. 

In what follows, we consider only stable circular orbits for which L 2 > 3c 2 r| 
implying, according to Eq. (17), that the radius of the orbit R, must be larger 
than 3r$. 



4 Energy quantization and Hawking tempera- 
ture 

As pointed out in Sec. 1, the quantization of energy is closely related to the 
quantization of the angular momentum of the orbiting test particle. This stems 
from the fact that the angular variable is a periodic function of time so that 
we may apply the prescription of Wilson [4] and Sommerfeld [5] 
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J, 



Ld(j) = nh 



(18) 



where denotes the action variable for one period of motion, and L is the 
canonical momentum conjugate to the angular variable (see Eq. (8)). 

Since L is a constant of motion, Eq. (18) yields the quantization condition 
for the angular momentum 

L = mL = nh (19) 

Equation (19) implies, via Eq. (17), a quantization of the radius R. For the 
smallest stable orbit, R = 3rg, the square root on the right hand side of Eq. 
(17) vanishes so that 



Rq 



f 2 
c 2 rs 



h 2 nl 
m 2 c 2 rs 



yielding 



where rc 



n 



3 2 r s 
r c 



(20) 



(21) 



^ is the Compton radius of the test particle. The next higher 
orbit corresponding to n\ = no + 1, has radius R\ given by Eq. (17) in the form 



r s I 



1 - 



3r| 



rhino + 1) 



(22) 



For the sake of simplicity, we assume that rs ^> rc- Eq. (21) then implies 
that no » 1 and (n + l) 2 ~ tIq(1 + — ) Expanding Eq. (22) in a small parameter 
2/n , we obtain 

i?l ~ Rq 1 + 

The true energy E = mE is obtained from Eq. (14), yielding 



no J 



(23) 



(24) 



From Eq. (16), we obtain 



R 2 



2R/r s - 3 

Substituting Eq. (25) into (24) gives the true energy as a function of R 

2 1 ~ rs/R 
E = mc — r 

(1 - lr s /R)i 
For Rq = 3rs, we obtain from this equation 



(25) 



(26) 
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= = .94mc 2 (27) 

3(0.5)3 v ' 

For given by Eq. (23), the true energy is given by 

*-*£tez(S P.) 

(1 + X — X^)2 

where x = (2/n )?. Expanding in small parameter x, we obtain from Eq. (28) 
to order x 2 

E 1 ^E {1 + ^-) (29) 
4n 

Then the excitation energy E\ — Eq = 8E\q is obtained from Eqs. (21), (27) 
and (28) as follows 

0.94 mc 2 ?ic 3 
^o-——- 0.068— (30) 

This formula shows remarkable resemblance to the expression for the tempera- 
ture of the Hawking thermal radiation [6] 

This similarity can be understood by examining the period of the angular 
variable 4>(t) for the particle encircling the black hole, and comparing it with 
the period of the Euclidean time te — it. For a circular orbit of radius R, the 
"orbital" period is given by 

2ttR 2ttR 2 . 
r = = -^—t (32) 

where we used Eq. (8) to write = ^f- = jfe. Using Eq. (25), we can express 
L/R in terms of R, and obtain from Eq. (32) 

2nR/2R \2/ 3 
r = — (--3) (l-- 2 rs/R) (33) 

where we used Eq. (A3) for i. Substituting the radius of the smallest stable 
orbit, R = 3rs, Eq. (33) yields 

T = 62 7T— (34) 

c 

Next, we consider the period of the Euclidean time te- Letting t e = it in 
Eq. (4), we obtain the Euclideanized metric. Following Ref. [10], we require 
this metric to be regular at the horizon. Thus, we define a small quantity p by 



r = r s (l + P 2 ) (35) 
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and expand the metric about p = 0. In this way, we obtain in the vicinity of 
the horizon 



ds 2 = 4r 



(36) 



The first two terms of Eq. (36) correspond to a metric dp 2 + p 2 d<p 2 in polar 
coordinates p, <p in flat plane. Since <f> has period 2tt, Eq. (36) implies that the 
period of te is 



Except for the dimensionless prefactor, the similarity to the orbital period (34) 
is apparent. 

The function 4>(te) also exhibits a period hp which is obtained by considering 
the partition function Z = T r (e~° H ) where (3 = (fc^T) -1 . For particle on a 
circle, Z can be expressed as a path integral with topological constraint making 
= and (j) — 27r indistinguishable. Z is then given as a sum of path integrals 
involving paths with various winding numbers Z [11] . These paths satisfy the 
condition 4>{te = h/3) = 4>{te = 0) + 2nl. In view of the topological constraint, 
path differing by 2irl are indistinguishable, implying that the period of 4>(te) is 
te = hp. Combining this result with Eq. (37), we obtain 



which agrees with Eq. (31). 

The fact that 5E and UbTh are quantities of the same order of magnitude 
suggests that interaction of photons with matter may play role in establishing 
the thermal equilibrium of the black body radiation [12]. For matter with quan- 
tized energy the transitions n — > n ± 1 correspond to absorption and emission of 
photons by matter. These effects should be particularly strong when 8E = hu) m 
= 2.8 ks Tjj, where ui m is the frequency at which the Planck formula for spec- 
tral distribution of black body radiation is maximum [12]. Using Eq. (31), we 
obtain huj m ss 0.1-p^r, which is in reasonable agreement with 8E\q given by Eq. 
(30). This implies that matter orbiting the black hole may be instrumental in 
establishing the thermal equilibrium of the radiation due to Hawking evapora- 
tion. There is another more direct way to demonstrate the relation between the 
dynamics of an orbiting particle and the Hawking radiation. This is the Unruh 
effect discussed in the Appendix. 

Now, we consider large orbits such that R 3> rg. 

From Eq. (17), we have 



Using the condition of mechanical stability, L 2 = MGR, the second term in 
the parenthesis can be written as 



47T 

te = — r s 



(37) 




(39) 
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2,c z ri 6r s 



< 1 (40) 



L 2 R 

Thus, to lowest order in rs/R, the radius of the orbit becomes 

if 2 

R^— (41) 

c rs 

To obtain the energy as a function of R, we start from Eq. (26). Expanding 
this equation in a small parameter rs/R, we obtain to lowest order in rs/R 

E~mc>(l-%L) (42) 
Introducing the radius R from Eq. (41), Eq. (42) yields 

The quantized energy E n follows from Eq. (43) by replacing L by quantized 
angular momentum L n = nh/m. In this way, we obtain 

E ^ mC \ l -S^) (44) 

The dependence of E n on the quantum number n is reminiscent of Bohr's quan- 
tization of the total energy of an atomic electron as expected for a Kepler 
problem. Using this result, we obtain the energy splitting 



SE n — E n+ i - E. 



mc 2 r% r 1 1 



n+l ~ 8f>2 



C 



' 2 (n+l) 2 



(45) 



For large values of n, the bracket can be replaced by 2/n 3 so that 



m 3 r|c 4 fiM 2 G 2 

where n has been eliminated using Eq. (19). Invoking the mechanical stability 
condition Z 2 = MGR, we obtain from Eq. (46) 

R2 H 

where the second equality stems from the relation v 2 = MG/R which is equiv- 
alent to the mechanical stability condition as seen by writing L = vR. Incin- 
dcntally, this result agrees with the Heisenberg uncertainty principle for energy 
and time. 

However, this Newtonian result becomes modified outside the galaxy core 
where v(R) tends to a constant value v . This has been discovered by measuring 
the Doppler shift in the 21 cm-line of neutral hydrogen in the galaxy disc [13]. 
This behavior can be attributed to the fact that there is a mass distribution 
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M(R) « 4? 



(48) 



where M(R) is the mass within the radius R. Hence, the corresponding mass 
density goes as p(r) oc r~ 2 indicating that there is a "hallo" of dark matter which 
stretches well beyond the visible part of the galaxy [13] 

5 Observational considerations 

In this section, we want to distinguish between "ordinary" black holes (OBHs) 
and primordial black holes (PBHs) . The OBHs are formed by collapse of massive 
stars (perhaps 10 — 10 10 times heavier than our Sun). The much lighter PBHs 
may have been formed by collapse of the overdense regions shortly after the Big 
Bang [8]. They are an important potential dark matter candidate [14, 15]. 

We first consider the supermassive black hole at the center of our own Milky 
Way galaxy. Its mass, obtained by carefull astrometric measurements [16], is 
M « 4 x 10 6 M Q where M = 1.989 x 10 33 g is the solar mass. The corresponding 
Schwarzschild radius is r$ ~ 1.2 x 10 12 cm. For the smallest stable orbit of radius 
R = 3r s , Eq. (30) yields SE W w 2.1 x 10~ 18 cV. Owing to the large mass of the 
black hole this quantity is extremely small. Clearly, the OBHs are not expected 
to produce a significant increase of SE W since their mass scale is of order 10 
Mq or more. 

On the other hand, the PBHs can be produced at a mass scale much smaller 
than Mq. For example, a black hole of mass equal to that of the Earth « 6 x 
10 27 g can be produced at time r w 10 _n sec after the Big Bang [8]. With the 
corresponding Schwarzschild radius r$ ~ 0.9cm, the formula (30) yields 8E W w 
2.8 x 10 - 6 cV. 

Let us assume that a test particle with mass of the hydrogen atom is orbiting 
this black hole at the radius 3r$- The corresponding Compton radius rc = ^ 
for this test particle is rc ~ 2 x 10 ~ 14 cm. The relevance of rc comes from the 
fact it gives the distance over which a relativistic particle of mass m is spread 
out about its position [9, 17]. For a test particle moving in the gravitational 
field of the black hole, the potential seen by the test particle becomes effectively 
averaged over a region of size rc ■ 

If rc <C r$ , the test particle behaves as a classical point mass moving in the 
potential given by Eq. (4). In this limit, the present semiclassical quantization 
approach is valid. 

On the other hand, as rc grows and approaches the value of r§, the present 
approach ceases to be valid and should be replaced by solving the complete 
quantum mechanical problem of Klein-Gordon scalar field in the space-time of 
the black hole [8]. 

The criterion of validity of the present approach, rs ^ rc, can be written 
as follows 



2mM > Mp 



(49) 
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where m is the mass of the test particle, M is the black hole mass, and Mp is 
the Planck mass M P = {hc/G)^ « 1.2 x 10 19 GeV/c 2 w 2.2 x 10 ~ 5 g. 

For the above considered case of hydrogen atom in the gravitational field 
of a black hole with mass equal to that of the Earth, the inequality (49) is 
well satisfied. However, a problem arises for PBHs that are lighter than about 
M rts 10 15 g. The importance of these black holes derives from the fact that 
they would shed all their mass by Hawking evaporation by now [18]. Inserting 
M « 10 15 g into Eq. (49), we see that the rest mass of the test particle should 
be much larger than about 2.4 x 10 ~ 25 g. This implies that our approach is 
marginally valid for m equal to that the hydrogen atom « 1.67 x 10 ~ 24 g, but 
fails for particles such as mesons and especially for neutrinos. For such particles, 
the energy quantization and its role in equilibration of the emitted black body 
radiation, discussed in the previous section, may not be realized. Obviously, as 
M goes well below 10 15 g, Eq. (49) imposes even stricter limits on the mas of 
the test particle. 

Now, let us consider the case of R ^> rs- The most accurate available 
information needed to evaluate the energy splitting from Eq. (47) comes from 
studies of the 21 cm - line of neutral hydrogen in our galaxy [19, 20]. From the 
wide spread of the data in the galactic latitude, Muller and Oort [20] deduce 
that the detected hydrogen cloud must be relatively close to the Earth. Hence, 
we use in Eq. (47) a value of R given by the sun-galactic center distance Rq w 
3 x 10 22 cm [7]. For v a , we use the sun-orbital velocity v « 2.2 x 10 7 cm/sec 
[21]. With these values, Eq. (47) yields an energy splitting SE n w 4.6 x 10 ~ 31 
eV. This result should be compared with the width of the hydrogen emission 
line [19] that has been measured to be about 80 kHz (which corresponds to 
an energy width of SE w 3 x 10~ 7 cV). Obviously, the energy quantization 
predicted by Eq. (47) is too small to be detected via studies of the 21 cm-line 
for hydrogen clouds orbiting the supermassive black hole at the center of our 
Milky Way galaxy. 

On the other hand, for the PBH with mass equal to that of the Earth, 
Eq.(30) yields an energy splitting SE W ~ 2.8 x 10 ~ 6 eV that is one order 
of magnitude larger than the width of the hydrogen emission line leading to 
possible observable effects. 



A Unruh effect for orbiting test particle 

According to Unruh [23] , an accelerated mass experiences itself to be imbedded 
in a hot gas of photons with a temperature Tjj given by 

k B T v = ^- (Al) 
2ttc 

where a is the proper acceleration. We now consider a test particle orbiting the 
black hole and calculate a for R given by the smallest stable radius R — 3rs- 
For a circular orbit, a is given by 
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a = a Q i (A2) 
where i — dt/dr is obtained from Eqs (7) and (26) in the form 

1 = c 2 (l-r s /R) = (l-lr s /R)i m 
The " ordinary" acceleration a is given by 

L 2 

Using Eq. (16), we can express the quantity L? in terms of the radius of the 
orbit as follows 

c 2 R 2 



Introducing this result into Eq. (A5), we have 

*o(R) = 7T~~ s (A6) 

-3) 

Eqs. (A2), (A3), and (A6) yield the proper acceleration for the orbiting test 
particle in the form 

a(R) = — {At) 

3„ 



2R 2 [i- irs/R 

In distinction from Ref. [23] where a(R) diverges as R — > r$, Eq. (A7) 
diverges as R — > Zrs/2. This is not surprising, since we are considering an 
orbiting test particle as compared with the "standing still" particle in Ref. [23]. 

1 2 

For R = 3r s , Eq. (A 7) yields a = Using this result in Eq. (Al) and 

invoking Eq. (38) we obtain 

k B T v = -(—^Mlk B T B (AS) 

This result and Eq. (30) suggest that the excitation energy SE W and the 
proper acceleration a(R) for R = 3rg are related quantities. To see this, we first 
express 8E\q in terms of the orbital period t as follows 



^ 10 = — -0.136^ (A9) 



where Eq. (34) has been used on the right hand side. This result is in good 
agreement with Eq. (30). Consequently, we may regard Eqs. (33) and (34) as 
suitable starting point for further discussion of proper acceleration. 
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From Eq. (A7), we obtain a ~ 0.157c 2 /r s . Using this result and Eq. (A9), 
we have 

a ~ 1.155^10 1 (-410) 
Substituting this result into the Unruh formula (Al), we obtain 

5E W ~ ^k B T v ~ 1.64fc B T ff (All) 
1.15 

where Eq. (A8) has been used to obtain the second equality. If we use Eq. (38) 
for k B T H , the right hand side of Eq. (All) becomes equal to ~ 0.13^ 

in good agreement with Eq. (30). This confirms the idea that the relation 
between SEi and the Hawking temperature T B discussed in Sec. 4 is actually 
the Unruh effect in disguise. 
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